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ELECTROMAGNETIC  BACKSCATTERING  FROM  A 
TRIANGULAR  DIELECTRIC  CYLINDER 
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IRE  NASU,  ul.  Proskury  12,  Kharkov  61085,  Ukraine.  E-mail:  lyusi@vil.com.ua 

INTRODUCTION 

Boundary  integral  equations  (lEs)  are  an  accurate  and  versatile  tool  in  the  electromagnetic 
scattering  by  smooth  cylinder  [1],  especially  when  it  is  combined  with  analytical 
preconditioning  [2],  An  original  approach  to  the  E-wave  scattering  by  PEC  polygons  was 
developed  in  [3]  to  study  a triangular  prism.  Here  the  boundary  of  the  scatterer  was  conformally 
mapped  on  a circular  cylinder  [4],  and  a log-singular  IE  was  treated  by  a simple  moment 
method.  We  combine  this  idea  with  the  method  of  analytical  regularization  (MAR)  and  study 
the  scattering  by  triangular  dielectric  cylinders.  Applications  of  this  analysis  are  expected  in  the 
design  of  microwave  and  optoelectronic  prism  sensors  and  couplers. 
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FORMULATION  AND  INTEGRAL  EQUATIONS 

The  geometry  of  the  problem  is  shown  in  Fig.l.  A plane 
electromagnetic  wave  //'=<? 
uniform  dielectric  cylinder  whose 
isosceles  triangle  of  the  base  2d  and  height  /?,  having 
relative  permittivity  e.  The  field  functions  have  to  satisfy 
the  Helmholtz  equation  in  free  space  and  in  dielectric,  for 
r g Q(r  eex  t Q)  and  re  Q , respectively.  On  the  contour 
L the  transmission  conditions  should  be  satisfied.  The 
scattered  field  must  also  satisfy  the  Sommerfeld  radiation 
condition  at  infinity'  and  condition  of  local  power  finitness. 
We  will  seek  the  solution  of  the  problem  in  the  form  of  a 
pair  of  single-layer  potentials  with  unknown  densities  ip  e C(d Q)  and  \\l  e C(d Q)  By  using 
the  conditions  on  L and  the  properties  of  the  normal  derivative  of  a single-layer  potential  when 
crossing  the  contour  of  integration  [1],  the  following  set  of  singular  lEs  is  obtained: 
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Fig.  1 : Shattering  geom etry  and  notations 
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To  facilitate  building  a fast  algorithm,  which  is  able  to  solve  the  scattering  by  cylinder  with  a 


contour  arbitrarily  close 
approximation  of  curve  L. 

, . d f N 

m=- 

M 


to  triangle,  we  shall  use  the  following  continuous  parametric 
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One  can  obtain  these  formulas  by  conformal  mapping  of  the  region  outside  a regular  polygon 
from  Z-  plane  (Z=x+iy)  to  the  (t,s)  plane,  where  0 < t < 2n  , and  by  using  the  transformation 

formula  first  derived  in  [5]:  dZ  / dt  = C(cos(3t  !2)\  ’,  where  C is  a complex  coefficient 
depending  on  the  size  and  orientation  of  the  polygon.  If  K-^> co  then  the  curve  given  by  (3)  tends 
to  a triangle,  however  by  truncating  the  series  in  (3)  at  the  W-th  terms  a smoothed  contour  is 
obtained  having  the  corner  radius  of  curvature  b inverse  proportional  to  K.  By  changing  the 
variables  in  (1)  to  ! and  t\  we  arrive  at  lEs  as  follows  (see  also  [2]): 
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| p{t')K e (t,t')dt'  - \q{t')K{t,  t')dt'  = f(t) 
0 0 


Here  the  following  notations  have  been  used: 


K(t,t')  = ~Ho)(kR(t,t')), 


Mz  (t,t ')  = — /(0~  — — x'(0 
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R(t,t')  = J(x(t)-x(t'))2  + (y>(/)  - .y(O)2  *s  the  distance  between  two  points  on  L.  Besides, 
in  the  right-hand  parts  we  have: 

/ (0  = exp  (ik(x{t)  cosy  + y(t)  siny  ))  (9) 

f (t)  = k(y'(t)  cos y — x'(^) siny )exp(/r(x(/) cosy  -t-j(r)siny))  (10) 


REGULARIZATION  AND  DISCRETIZATION 

As  a curve  given  by  (3)  is  smooth,  we  can  solve  (5)  by  projecting  them  onto  the  set  of  global 

{init  l+co 

e oo  • Here  the  key  step  is  to  split  the  IE  kernels  by  adding  and 

subtracting  the  terms  corresponding  to  canonical  scatterer,  i.e.,  a circular  cylinder  [2].  Thanks  to 
the  fact  that  the  singularities  of  actual  kernels  are  the  same  as  of  canonical  ones,  and  that 
angular  exponents  form  a complete  set  of  orthogonal  eigenfunctions  of  the  canonical  operators, 
the  resultant  infinite-matrix  equation  is  of  the  Fredholm  second  kind.  It  has  a (2x2)  block 
structure  generated  by  (5)  in  obvious  manner.  Therefore  this  specialized  Galerkin  projection 
procedure  can  be  regarded  as  an  analytical  preconditioner,  which  leads  to  an  always-stable 
numerical  solution  whose  accuracy  is  easily  understood  and  controlled.  Namely,  if  the 
intermediate  computations  have  been  done  with  superior  accuracy,  then  the  error  in  matrix 
inversion  is  controlled  by  the  size  N of  each  block,  so  that  large  enough  N guarantees  a uniform 
accuracy.  An  important  point  of  the  MAR-algorithm  efficiency  is  fast  computation  of  the  matrix 
elements,  which  are  the  double  Fourier  transform  coefficients  of  the  twice-continuous  functions 
formed  by  the  differences  of  actual  and  canonical  kernels  of  IEs.  Here  various  versions  of  FFT 
and  the  like  numerical  algorithms  can  be  successfully  implemented.  We  shall  demonstrate  the 
features  of  the  developed  method  by  presenting  the  dependences  of  the  matrix  condition  number 
and  inversion  error  as  a function  of  the  matrix  block  size  N for  various  triangle  sizes,  dielectric 
constants,  peak  curvatures  of  smoothing  parameterization,  and  frequencies. 
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